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We consider the level 1 solution of quantum Knizhnik-Zamolodchikov equation with re- 
flecting boundary conditions which is relevant to the Temperley-Lieb model of loops on 
a strip. By use of integral formulae we prove conjectures relating it to the weighted enu- 
meration of Cyclically Symmetric Transpose Complement Plane Partitions and related 
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1. Introduction 



Since tlie papers |ll|Jl], tliere lias been a great deal of work on the combinatorial 
interpretation of quantum integrable models at special points of their parameter space. 
The original observation is that the numbers of Alternating Sign Matrices (ASM) and Plane 
Partitions (PP) in various symmetry classes appear naturally in the ground state entries 
of the Temperley-Lieb 0(r = 1) model of (non-crossings) loops with various boundary 
conditions (and related models) . The appearance of ASM numbers was developed further 
and to some extent explained by the Razumov-Stroganov conjecture ^ and variants 
interpreting each ground state entry as a number of certain subsets of ASM. The role 
of plane partitions remained more obscure until the recent work [^|^ which showed that 
the enumeration of symmetry classes of PP also occurs naturally on condition that one 
consider a slightly more general problem, namely the quantum Knizhnik-Zamolodchikov 
equation (gKZ), first introduced in this context in and in which the parameter r is 
now free. This provided a (conjectural) bridge between enumerations of symmetry classes 
of ASM and PP, which is a fascinating topic of enumerative combinatorics in itself. 

The present work is concerned more specifically with the case of the Temperley- 
Lieb loop model (and its qKZ generalization) defined on a strip with refiecting boundary 
conditions (the case of periodic boundary conditions was treated similarly in |^]). The 
corresponding ASM were discovered in [^|10[' they are Vertically Symmetric Alternating 
Sign Matrices (VSASM) of size (2n-f 1) x {2n + 1) in even strip size = 2n, and modified 
VSASM of size (2n + 1) x (2?! + 3) in odd strip size = 2n + 1. As to the PP, they were 
discussed in [0: they are Cyclically Symmetric Transpose Complement Plane Partitions 
(CSTCPP) in odd strip size, and certain modified CSTCPP (referred to as CSTCPP^ 
in the following) in even strip size. The conjectures of concerning the r-enumeration 
of these plane partitions are the main subject of this work. In Sect. 2 we shall review the 
basics of integrable loop models based on the Temperley-Lieb algebra; in Sect. 3 we shall 
discuss the related qKZ equation, and review the conjectures of 0; in Sect. 4 we introduce 
the main technical tool, that is certain explicit integrals solving qKZ; and finally in Sect. 5 
and 6 we prove the conjectures of 0, considering separately even and odd cases. 
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2. Loop model with reflecting boundary conditions and link patterns 

2.1. Dense Loop model on a strip 

We consider the version with reflecting boundarieslll of the inhomogeneous 0{1) non- 



crossing loop model [|T2[. The model is defined on a semi-infinite strip of width (even 
or odd) of square lattice, with centers of the lower edges labelled 1, 2, A^. On each face 
of this domain of the square lattice, we draw at random, say with respective probabilities 
1 —ti,ti in the z-th column (at the vertical of the point labelled i) one of the two following 
configurations 



f2.1) 



The strip is moreover supplemented with the following pattern of fixed configurations of 
loops on the (left and right) boundaries: 
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With probability 1, a configuration will lead to a pairing of the points 1,2, ...,A^ 
according to their connection via the paths (except for one point if is odd which is 
connected to the infinity along the strip). Such a pattern of connections is called a link 
pattern, and an individual pairing is called an arch. The set of link patterns on points 
is denoted by LP^, and has cardinality (2n)!/(n!(n + 1)!) for A^ = 2?! or A^ = 2n — 1. 
Each link pattern of odd size A^ = 2?! — 1 may indeed be viewed as a link pattern of size 
2n but with the point 2n sent to infinity on the strip: this provides a natural bijection 
between LP2n-i and LP2n- A link pattern n G LP^ may also be viewed as a permutation 



^ These boundary conditions are sometimes called "open" , in reference to the equivalent open 
XXZ spin chain, or "closed" , due to the way the loops close at the boundaries of the strip. 
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Fig. 1: A sample configuration of the Dense Loop model on a strip of width 
= 6 (left). We have indicated the corresponding open non-crossing link 
pattern of connection of the points 1, 2, 3, 4, 5, 6 (right). 

TV E Sn with only cycles of length 2 (except one cycle of length 1 for odd), and we 
shall use the notation n{i) = j to express that points i and j are connected by an arch. 
For a pair (i, j) such that j = n{i) and i < j, we will call i the opening and j the closing 
of the arch connecting i and j. An example of loop configuration together with its link 
pattern are depicted in Fig. 1. We use a standard pictorial representation for link patterns 
in the form of configurations of non-intersecting arches connecting regularly spaced points 
on a line, within the upper-half plane it defines. For odd A^, the unmatched point may 
be represented as connected to infinity in the upper-half plane via a vertical half-line. 
We moreover attach a weight t = —{q + q~^) = 1 to each loop (hence the denomination 
0(r = 1) model, q = -e^^/^). We may then compute the probability Prob(7r) for a given 
randomly generated configuration of the loop model on the strip to be connecting the 
boundary points according to a given link pattern tt. 

In Ref. [|12|, the model was solved by means of a transfer matrix technique, using 



solutions of the boundary Yang-Baxter equation [13| \T^] that parametrize the inhomoge- 
neous probabilities ti via integrable Boltzmann weights, coded by a standard trigonometric 
i?-matrix. Using the integrability of the system, and following the philosophy of [^, the 
suitably renormalized vector of probabilities \1/ = {\l/7r}7r£^p^ was shown to satisfy the 
quantum Knizhnik-Zamolodchikov equation with refiecting boundaries for q = — e*'^/'^, in 
the link pattern basis. 

In the following, we will consider the more general case of generic g, r, which does 



not have stricto sensu an interpretation in terms of lattice loop model |16 
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2.2. R matrix 

The i?-matrix of the model is an operator acting on hnk patterns of LP^: 

qw-q ^z qw-q ^z 

qz — q~^w z — w 

= —I H — e^ 

qw — q ^z qw — q ^z 

where z and w are complex numbers attached to the points i and i -\- 1 and where e^, 
i = 1, N — 1 are the generators of the Temperley-Lieb algebra TLj^{t), subject to the 
relations 

ef=Tei, [ei,ej] = if |z - j| > 1, €^€^±16^ = 6^ (2.4) 
with the parametrization 

r = -g-g-i (2.5) 




/■ i+1 



Fig. 2: Action of the Temperley-Lieb generators on link patterns. 

In ( p.3| ), we have depicted the Temperley-Lieb generators as tilted squares with edge 
centers connected by pairs. The corresponding action on link patterns should be under- 
stood as follows (see Fig. 2): assume the points + l are connected to say the points j, k 
in a link pattern tt. Then, unless j = i + 1 and k = i, the link pattern n' = e^Tr is identical 
to TT except that i is now connected to z + 1, and j to k. If j = i + 1 and k = i, the points 
i,i + 1 are connected to each other in tt, and e^Tr = ttv. The latter is a direct consequence 
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of the projector condition e? = re^, as any link pattern with i connected to z + 1 hes in 
the image of e^. 

The above i?-matrix satisfies the Yang-Baxter equation and the unitarity relation 



Ri^i+i{z, w)Ri+i^i+2{z, x)Ri^i+i{w, x) = Ri+i^i+2{.w, x)Ri^i+i{z, x)Ri+i^i+2{z, w) 
Ri^i+i{z, w)Ri^i+i{w, z) = I 



as consequencs of the Temperley-Lieb algebra relations (27i). 



(2.6) 



2. 3. Link patterns, dyck paths, and containment order 




(a) 



(b) 




(c) 



Fig. 3: Dyck path (b) associated to a link pattern (a). The former is obtained 
as the discrete path on the non-negative integer line: 0, 1, 2, 1, 2, 3, 2, 1, 0, 1, 0. 
We have represented in (c) the box decomposition of the Dyck path. 



Before turning to the gKZ equation, we wish to emphasize a number of useful prop- 
erties satisfied by the link patterns, and the action of the Temperley-Lieb generators. An 
alternative pictorial representation of link patterns is via Dyck paths, namely paths from 
and to the origin on the non-negative integer line, with steps of ±1 only. The bijection 
between link patterns and Dyck paths is illustrated on Fig. 3 for = 10. For even N = 2n, 
we construct the Dyck path by visiting the points of the link pattern from 1 to N, starting 
from the origin of the non-negative integer half-line, and with the following rule: if an 
arch opens (resp. closes) at z, then the path goes up (resp. down) one step between time 
i — 1 and i. The path is guaranteed to come back to the origin at time 2?! as there are as 
many openings as closings of arches, and moreover stays in the non-negative half-line as 
all arches must first open before closing. In the case of odd N = 2n — 1, one arch exactly 
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has an opening and no closing point (it is connected to infinity) , hence the path ends up at 
the point 1 on the integer half-hne. It can be completed into a path of length 2n by a final 
step to the origin, thus expressing on Dyck paths the abovementioned bijection between 
LP2n-i and LP2n- The Dyck path is represented in the plane as the (broken-line) graph 
of the function (t, h{t)), t = 0, 1, iV. 

Dyck paths allow to endow the set of link patterns with a natural "containment" order, 
namely tt < tt' iff the Dyck path of tt contains strictly that of n' . This notion is made 
explicit by introducing the "box decomposition" of any given Dyck path (see Fig. 3 (c)), 
namely the decomposition of the region of the plane delimited by the path and a broken 

fine (0, 0) (1, 1) ^ (2, 0) ^ ^ (iV - 1, 1) ^ {N, 0) if = 2n, without the last step 

for A'' = 2n — 1, by use of squares of edge \/2 tilted by 45°. Then a Dyck path 5 contains 
strictly another 5' iff S is obtained from 6' by addition of at least one box. A box addition at 
position m consists simply in replacing a portion of path (m— 1, h) — * (m, h—1) (m+1, h) 
that visits successively the points h, h — 1, h of the integer half-line at times m — 1,171, m+l, 
by the portion (m — l,h) — > (m, h + 1) ^ (m -|- l,h), thus adding a box with center at 
coordinates {m,h). This may also be described as transforming a local minimum into a 
local maximum at position m on the path. The "smallest" link pattern (whose Dyck path 
contains all others) is the pattern ttq with links 7^o(^) = 2n+l — i, i = 1, 2, n for A?" = 2n, 
and i = 2, ...,n, for N = 2n — 1, while 7ro(l) = 1. It corresponds to the farthest excursion, 
reaching point n on the integer half-line. The "largest" link pattern (whose Dyck path 
is contained in all others) TTmax has 7rmax(^) = ■i + 1 for z = 1, 3, .., 2n — 1 when N = 2n 
and i = 1,3, 2n — 3 when N = 2n — 1, while 7rinax(2n — 1) = 2n — 1. It corresponds 
to the shortest range excursion, alternating betwen the origin and point 1 on the integer 
half- line. So we have tiq < it < TTmax for all tt G LPn such that tt 7^ ttq, TTmax- Finally, we 
shall denote by P{tt) the total number of boxes in the box decomposition of the Dyck path 
associated to tt. We have for instance /3(tto) = n{n — l)/2 and /9(TTmax) = 0. 

The action of on link patterns may be easily translated into the language of boxes 
on Dyck paths. The action of may indeed be viewed as a box addition at position i on 
the corresponding Dyck paths. Then 3 situations may occur (Fig. 4): 

(i) The path has a minimum at point i: the box addition transforms it into a maximum. 

(ii) The path has a maximum at point i: the box-added path is unchanged, but picks up 
a factor of r. 

(iii) The path has a slope at i, namely a succession of two up or two down steps: the box 
addition actually destroys the two rows of boxes at its height and immediately below 
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Fig. 4: Box addition at position i on Dyck paths corresponding to the action 
of Ci. We have depicted three generic situations for the addition: (i) at a local 
minimum (ii) at a local maximum (iii) at a slope. Both (ii) and (iii) lead to 
a Dyck path contained by the original one, while (i) produces a Dyck path 
containing it, with exactly one additional box. 

until the other side of the path is reached. The net result may be interpreted as an 
avalanche, in which the mountain shape between the point of impact and the other 
side falls down by two units. 

This allows to see that among all possible actions of on a link pattern tt, only one 
leads to a "larger" Dyck path (containing tt): Ciiv = n' < n, namely in the situation (z), 
while all other situations lead to n < n' = Cin. This observation will be used below. 

The interpretation of the action of on Dyck paths was used in 0] to rephrase the 
homogeneous loop model as the stochastic model of a growing interface. 

3. The qK.Z equation for reflecting boundary condition 

3.1. The equation 

The reflecting boundary qKZ equation consists of the following system of equations 
for a vector \E' which depends polynomially on the variables zi, . . . ,zn (and q~^): 

Ri{zi+i, Zi)'^Nizi, ...,Zi, Zi+i, ...,zn) = ^'Ar(^i, • ■ ■ , Zi+i, Zi,..., Zn) (3.1a) 
cn{zn)'^n{zi, . . . , Zn) = ^'Ar(^i, . . . , zl-1, s/zn) (3.16) 
ci{zi)'^n{zi, . . . , Zn) = ^'Ar(l/Zl, ^2, • • • , Zn) (3.1c) 
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Here ci and cn are scalar functions to be specified later, and s = ^^('=+2) is a parameter 
which determines the "level" k of the equation: here we consider the so-called level 1 case, 
namely with s = . 

One can think of Eqs. ( |3.1| ) as an analogue of the quantum Knizhnik-Zamolodchikov 
equation (gKZ) in the form introduced by Smirnov (see also fl^), in which one re- 
places the periodic boundary conditions, implicit in the usual qKZ, with reflecting bound- 
1l4| . More precisely, Eq. ( p.la|) is the exchange relation corresponding to the bulk. 



aries 



independent of boundary conditions, whereas Eqs. (|3.16, c|) implement the reflections at 
the two boundaries. 

In [^, it was remarked that solving these equations for even size N = 2n automat- 
ically provides a solution for odd size — 1 by taking the last parameter zn io zero 
(or equivalently to infinity). We therefore discuss in detail the case of even size now, 
postponing to Sect. 6 the discussion of the odd case. 

3.2. Minimal polynomial solution 



In [|T2|, it was claimed that the system of equations ( |3.1|) possesses a polynomial 
solution of minimal degree 3r2,(n — 1) which is unique up to multiplication by a scalar. To 
actually solve the equations ( |3.1|) , one first remarks that the — 1 equations ( p.la|) from 
a triangular system with respect to the containment order of Dyck paths introduced in 
Sect. 2.3. Indeed, when written in components, this equation reads: 

^-^^^^-^^{n-l)^^{zu...,Z2n)= Y] *-'(^l,---,^2n) (3.2) 
Zi+1 - Zi f—f 

ein'=TT 

where acts on functions of the 2;'s by interchanging Zi and Zi-^-i. Now consider the sum 
on the r.h.s.: it extends over the proper inverse images of n under e^. Picking tt in the 
image of (i.e. with an arch connecting points i and z -|- 1, as explained above), its inverse 
images tt' under all have dyck paths containing that of n (i.e. n' < n) except one, 
say TT*, corresponding to the Dyck path of n with the box at position i removed, hence 
with n < n*. Hence Eq. (|3.2| ) allows to express \E'7r* in terms of only \1/q,, with a < n*. 
The solution is therefore uniquely fixed by specifying the component corresponding to the 
smallest link pattern ttq defined above, whose Dyck path that contains all others. The 



latter is entirely fixed by the degree condition and factorization properties deduced from 
the qKZ system; the result is: 



= n (^^i^^ - ZiZj) Yl (Q^i-Q ^Zj){q'^ - ZiZj) (3.3) 

l<i<j<n n+l<i<j<2n 

This fixes the functions ci and cat in Eqs. ( p. 16, c| ) to be ci(x) = l/x'^"^"^ and cn{x) = 

It is also a simple exercise to prove that the solution enjoys a reflection invariance 
property, inherited from the inversion (unitarity) relation satisfied by the i?- matrix ( |2.6|) , 
and clearly satisfied by the fundamental component ^E'ttq: 




where p{n) is the mirror image of tt w.r.t. a vertical axis. 

At the special value q = — exp(z7r/3), which is a "degenerate" case of the qKZ system 
where the shift s becomes unity, the vector \E' can be interpreted as the ground state 
eigenvector of the loop model described in Sect. 2. For generic q, no such direct physical 
interpretation of \1/ exists; however, it still retains remarkable combinatorial properties that 
we shall describe now. 



3. 3. Conjectures 

In Ref. 1^, a number of conjectures were made on the homogeneous, generic q limit 
of the components of The dependence on q was always expressed through the quantity 
r = —q — q~^, and the entries of \E' were observed to be polynomials of r with non-negative 
integer coefficients. The main conjectures concerned an identification of the sum of com- 
ponents of \E' with generating polynomials for weighted rhombus tilings of finite domains 
of the triangular lattice, involving the catalytic variable r. These were interpreted in odd 
and even sizes as r-enumerations of respectively Cyclically Symmetric Transpose Comple- 
mentary Plane Partitions (CSTCPP), or equivalently of rhombus tilings of a hexagon with 
cyclic and transpose-complementary symmetries (odd size) , and in the even case as modi- 
fied CSTCPP (CSTCPP^) corresponding to rhombus tilings of a hexagon with a central 
triangular hole and with cyclic and transpose-complementary symmetries. These conjec- 
tures were extensions of earlier conjectures concerning only the particular value r = 1 
{q cubic root of unity), due to Razumov and Stroganov 0], and proved in [jl2|, involving 
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respectively the total number of CSTCPP in odd size, and that of Vertically Symmetric 
Alternating Sign Matrices (VSASM) in even size. 

Apart from the main sum rule conjectures, a number of other conjectures were made 
in ^ on the entries of \1/ as polynomials of r, concerning degree and valuation, and also 
conjecturally relating the small r behavior to the numbers of Totally Symmetric Self- 
Complementary Plane Partitions (TSSCPP) with fixed features, namely, once expressed 
as Non-Intersecting Lattice Paths (NILP), with fixed termination points of the paths. 

The purpose of the present paper is to prove these various conjectures, by means of 
multiple integral expressions for the homogeneous solution to the gKZ equation ^'(r). 

4. Integral expressions for solutions of level 1 qK.Z 

4-1- The intermediate basis 

The method introduced in in order to obtain integral representations of was to 
exhibit a different basis than that of link patterns in which the integral expressions for the 
components are relatively simple. We shall use the same basis here. Note that this section 
is "boundary conditions-independent" and its results are equally valid for say periodic 
boundary conditions. 

The elements of the basis we consider are indexed by strictly increasing integer se- 
quences of the form a = {ai, 02, . . . , a^}, where 1 < < 2z — 1. We denote by the 
set of such sequences. The sequences in On are in one-to-one correspondence with the link 
patterns in LP2n in two different, inequivalent ways. One may indeed associate to each 
TV E LP2n the sequence ai(7r), i = 1, 2, . . . , n, of the positions (counted from left to right, 
and taking values in {1, 2, . . . , 2n — 1}) of the openings of arches in n. Similarly, we may 
associate to tt the sequence 6i(7r) = ai{p{Tv)), z = 1, 2, . . . , n recording the closing positions 
of the arches of tt, counted from right to left, or equivalently the opening positions of the 
arches in the reflection p{n). 

In 0, we have constructed the change of basis from the link pattern basis to the 
"arch opening" basis, namely expressed the solutions \E'7r5 ^ LP2n in terms of multiple 
residue integrals \E'ai,...,a„, with {ai, 02, . . . , a^} G On- More precisely, Ref. expresses 
any integral \l/a = ^ai,...,a„ for weakly increasing sequences of a's as well, in terms of the 
solution \l/7r in the link pattern basis, via the linear transformation: 

^^izu...,Z2n) = Yl '^a,7r(r)*,r(-2l,...,22n) (4.1) 
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with polynomial coefficients Ca,7r(T) expressed as follows. Let Uk = {q'^^^ — — 
q~^). For k > 0, Uk is the Chebyshev polynomial of the second kind associated to — r, 
defined recursively by Uk+i = —rllk — Uk-i with Uq = 1 and U-i = 0. A given link 
pattern tt G LP2n may alternatively be thought of as a permutation of {1, 2, . . . , 2n} with 
only cycles of length 2. The arches forming tt may therefore be described by the pairs 
(z, 7r(z)), 1 < z < 2n, such that i < 7r(i), in which case i are the positions of the openings 
and 7r(z) of the closings of the arches in tt. Then we have the following formula: 

2n-l 

C^Ar)= n ^M(a,.) (4.2) 

i=l 
i<7r(i) 

where 

//(a, i) = card < aj < 7r(z)} H ^ (4-3) 

In 0, it was mentioned that the entries (^3) may easily be computed by iteratively 
removing the "little arches" of n (with n{i) = i + 1), and replacing them with a factor 
Um-i, where m is the total number of a's lying under that arch (namely such that aj = i). 
The new link pattern thus obtained has one less arch, and its a's are relabeled accordingly, 
while m — 1 extra a's are placed in position i — 1. The algorithm is iterated until the link 
pattern becomes empty. 

If we moreover restrict the set of a's to 0„, we get a true change of basis from 
link patterns to arch openings, in which C(r) is a square invertible Cn x Cn matrix with 
polynomial coefficients. The matrix C(r) indeed enjoys the following properties. Let us 
first use the bijection between LP2n and On to write Ca,7r(T) = CQ,7r(r), where a G LP2n is 
uniquely determined by its arch opening positions ai, a2, . . . a^, counted from left to right. 
Let us moreover order the link patterns, say by increasing lexicographic order on the sets 
of positions of their arch openings. Then we have the property 

(P) C(r) is a lower triangular matrix, with entries 1 on the diagonal, and polynomials of 

r with integer coefficients elsewhere, and the same holds for C~^(r). 
Property (P) is easily derived as follows. First, it is clear that the diagonal terms Ca(7r),7r = 
1, as /i(a(7r), i) = for all arch openings i of tt: indeed, the a'a being the arch openings of 
TT, the set {aj\i < aj < vr(z)} has one a per arch enclosed by the arch i 7r(z), henceforth 
a total of (7r(z) — i + l)/2. Consequently, all the indices of the Chebyshev polynomials 
contributing to (|4.2| ) vanish, and as C/q = 1? the result follows. The triangularity is best 
understood by following the abovementioned algorithm for constructing C. Indeed, at any 
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step in the algorithm, the only cause for the matrix element to vanish is that one has 
no a's under the little arch considered, as in this case one would get a factor U-i = 0. 
The matrix element Ca^-K can only be non-zero if the arch openings of a occupy positions 
lexicographically larger that those of tt. Indeed, by contradiction, assume the structure of 
arch openings in a and tt are the same up to say a position i where an arch opens in a 
while an arch closes in tt. This means that the total number of arch openings for positions 
j > i is strictly larger in tc than in a. Consequently, applying the above algorithm to the 
arches of tt opening at positions > z, we see that at least one little arch in the process will 
have no arch opening of a below it, thus receiving a weight U-i = 0, and therefore the 
corresponding matrix element C^^n vanishes. 

Finally, one can rewrite the qKZ equation itself using the linear combinations defined 
by Eq. ([4.1|) . Note here that we are forced to use not only the components corresponding to 
our basis On oi increasing sequences, but also those corresponding to any non-decreasing 
sequence. In principle all of them can be reexpressed as linear combinations of increasing 
sequences only, but it is preferrable to avoid having to write these linear dependence 
relations explicitly. 

All that is needed is the action of the on the \l/a- We have the following 
Theorem 1: For any non-decreasing sequence ai, . . . , such that the number i occurs 
exactly k times, k > 0, we have the formula: 



,2,. ..,2,. ..,(3,^ 

fc-1 



(ei\E')ai,...,i,...,i,...,a„ —Uk-lUk-4:'^ai,...,an — C^fc- 1 t^A;-3 (^ai , . . 1, 
fc 

+ ^E'ai,...,v^,i+l,...,a„) + Uk-lUk-2'^a^,...,i-l,i^^,i+l,...,a„ (4.4) 
fc-1 fc-2 

(where for k = the r.h.s. is zero). 

Proof: expand the l.h.s. in the basis of link patterns by using Eq. We find: 



(e2^')a= C^M^)'^^ + ^ Yl C'a,^^- (4-5) 

7r;7r(i):^i+l 7r:7i"(i)=i+l 

where we have distinguished among the entries of its diagonal entries, equal to r, and 
its non-diagonal entries, equal to 1. 

The \E'7r must be regarded here as independent objects, so that we must now check 
Eq. ( |4.4| ) for each link pattern tt. This will be performed by a case by case analysis of the 
situation around the sites + 1). Each time only the coefficients involving the arches 
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starting or ending at i, i + 1 differ from term to term in the equation, so that we can ignore 
the remaining factors. The proof wiU be explained pictorially using the same conventions 
as in appendix 1 of that is by drawing the coefficient Ca,7r as the usual (local) depiction 
of the link pattern tt decorated by placing between sites i and i + 1 (inside a circle) the 
total number k of a's such that aj — i. 
There are 4 cases: 

(i) If i is an opening and z + 1 a closing of tt, then tt has a little arch (i, i + 1): n{i) = i + 1. 
In this case the equality reduces pictorially to 

r /^~\ =Uk-iUk-4 /~\ -Uk-iUk-3 

■ • • * * • ■ • • ■ • ' 1 • • • • • • ^B- 



/ i+l i i+1 i i+1 

-Uk-iUk-3 +Uk-iUk-2 

! • . \ • ' • • • • ' • ! • . \ • • • ■ 

I i+i I i+i 

or explicitly 

rUk-i = Uk-iUk-A X Uk-i - Uk-iUk-3 X 2Uk-2 + Uk-iUk-2 x Uk-3 (4.6) 

which is easily checked by noting that Uk-iUk-4 — r = Uk-2Uk-2,- 
In all other cases there is no little arch (i, i + 1). 
(ii) If both i and i + 1 are openings, call p the total "weight" under the arch leaving i + 1, 
that is p = card{£|z + 1 < < 7r{i + 1)} — T^Ct+i) i ^ ^^^^ ^ ^^iq remaining weight 
under the bigger arch starting from z, excluding what is under the smaller arch and 
the weight k under the segment [i,i + 1), in order to make the pictorial description 
simpler: q = card{£|7r(i + 1) < ae < 7r(z)} — ""(^)~^(^+^)~^ _ Then the identity to prove 
is: 




Uk-lUq-l —Uk-lUk-4 X Up-lUk+p+q-2 — Uk-lUk-3 X Up-lUk+p+q-3 

— Uk-lUk-3 X UpUk+p+q-2 + Uk-lUk-2 X UpUk+p+q-3 (4.7) 

which is again a routine check. 
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(ii') The case where i and i + 1 are both closings is treated analogously. 

(iii) Finally, if z is a closing and z + 1 an opening, call p the weight under the arch (7r(z), i) 

defined as before, and q the weight under the arch (i + 1, 7r(z + 1)). Similary the proof 

of the identity 




i i+l 



Uk-lUk-A 
Uk-lUk-3 



Uk-iUk- 



i i+l 



i 1+1 




Uk-lUk+p+q-2 =Uk-lUk-4 X Up-lUq-l — Uk-lUk-3 X UpUg-l 

- Uk-iUk-3 X Up-lUg + Uk-iUk-2 X UpUq 
is left to the reader. This completes the proof of the theorem 1. 



(4. 



4-2. Integral solution of qKZ equation: general principle 

The idea to use integral representations for solutions of the qKZ equation is not new 
and there is a vast literature on the subject (cf the references in Sect. 11.2 of |]l^). We 
consider here a very specific type of level 1 solutions, for which one expects a much simpler 
formula than generically. In the present context, this idea was used in in the case 
of the qKZ equation with the usual periodic boundary conditions. We now describe the 
procedure in a slightly more general (boundary conditions-independent) setting. 

The idea is to define for any non-decreasing sequence (ai, . . . , a„) the following quan- 
tity: 



'^a,,...,a„{zi,...,ZN) 

F{zi, . ..,zn;wi 



q Zj 



n 



dwi 



2m 



n (^^^ 

l<i<j<N " " m=l 

Ul<e<m<ni'^n^ - Wl){qWi 



Wn): 



q ^Wr 



lV^=iUZiiwe-z,)YlZae+i(^^i-~1 ^ 



(4.9) 



where F is any rational function that is symmetric in all Zi and symmetric in all wi. The 
contours of integration encircle the Zi but not q~'^Zi, nor any poles of F. 

We wish to prove that ( [4. 91 ) solves the exchange relation ( |3.2D of the qKZ equation, in 
the operatorial form = (e^ — r)\I/ where the action of is given by Eq. ( [1.4| ), and we 
have introduced he divided difference operator ti = {qzi — q~^ Zi+i)di. As before, we have to 
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consider the cases where exactly k a's take the value i, say a^, = am+i = ■ ■ ■ = am+k-i = h 
while ttm-i < i or m = 1 and am+k >i or m + k = N. Starting from the expression ( [4.4| ), 
and subtracting r\E' on both sides, we wish to express (e^ — r)\E' as the action on ^' of the 
divided difference operator ti. We are left with proving the following 
Theorem 2: The function ( \i.9{) solves the exchange relation of the qKZ equation, namely 
it satisfies: 



fc 



{ti^)ai,...,arn-i,i^^^jji^,arri+k,---,a.n ~ Uk — 2Uk—3,^ a^,. 

k 

~ Uk-lUk—s{^ai,...,am-i,i — l,h^^,am+k,---,an ^ai ,.. .,a^_i ,v^.i+l,"m+fe ) 

fe-i fe-l 

+ Uk — lUk — 2^ai,...,a^-i,i — l,i,...,i,i+l,am+k,---,<^n 



(4.10) 



Proof: Two important remarks are in order. Firstly, the operator ti acts only on the pieces 
of \1/ that are non-symmetric in [zi, ^i+i). When acting with ti on ( [4.9| ), we may restrict 
our attention to the non-symmetric part of the integrand. Secondly, we note that for any 
function S{ui, . . . , Uk) satisfying the following vanishing antisymmetrizer property that 

5^(-l)^^K(i),...,M,(fc)) = (4.11) 

o-eSk 

dukS{ui,. . .,Uk) W {um-ui)=0 (4.12) 

l<e<m<k 

Indeed, for any given permutation a G S"/;;, we may perform the change of variables Um = 
v^^rn) for m = 1,2,..., A;, resulting in h = f dvS{v'^)A{v'^) = f dv{-iy S{v'')A{v), 
where we have used the antisymmetry of the Vandermonde determinant, and finally sum- 
ming over all permutations yields ( ^4.12| ), thanks to ( |4.11j ). We also note that if 5" satisfies 
( fl.lll) , any symmetric function of the w's multiplied by 5" will also satisfy it. 

As a consequence of the two above remarks, to prove the identity ( [4.10|) , it is sufficient 
to prove a weaker statement on the part P of the integrand of \1/ that is non-symmetric 
in {zi, Zi+i) and also non-symmetric in {wrmWm+i,---,Wm+k-i), after factoring out a 



then the multiple integral 
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Vandermonde determinant of the w's. Rewriting Ui = Wm+i-i for z = 1, 2, . . . , /c, P reads 
simply 

Pk = {qz,-q — (4.13) 

Now introducing the quantity 

S{ui, ...,Uk)= UPk - (Uk-2Uk-3 + Uk-iUk-sifi +9k) - Uk-iUk-2fi9k^Pk (4.14) 

where we use the notations 

f Ui-Zi qui-q-^Zi+i 

J I = n — 9t- = (4.15) 

qut - q ^Zi ue- Zi+i 



for £ = 1, 2, . . . , /c, we are simply left with the task of proving that (|4.11| ) is indeed satisfied, 
and ( [4.101 ) will follow from the above considerations. 

For simplicity, we will work with the quantity T = Sx 11^=1 (^^ — Zi){qui — 

q~^Zi+i), proportional to S of ( |4.14D by a factor symmetric in the u's, henceforth proving 
A{S) = amounts to proving A{T) = 0. Explicitly, we have 



a(t{ui, . . . , Uk)j = A(^Aq{u)^ i (qzi+i - q ^Zi) Y\ ftgi ~ (QZi - q ^Zi+i) 

\ i=i 

- {Zi+l - Zi)A(^Aq{u){Uk-2Uk-3 - Uk-lUk-sifl +gk) + Uk-lUk-2figk] 



(4.16) 



where we have written for short the Vandermonde as Ag(u) = Y[i<£<m<k(^'^^ ~ ?~^^m) 
and noticed that only this piece of P is non-symmetric in the w's. 

We need the following lemma, expressing the antisymmetrization of the g- Vandermonde: 

^(a,(u)) = {-l)'^('^-'y^U^U2 ■ ■ ■ Uk-i A(u) (4.17) 

Proof. We proceed by induction. For /c = 2, an explicit computation leads to A{qui — 
q~^U2) = {q + q~^)iui - U2) = -t/iA(u). Assume ( |4.17| ) holds for k ^ k - 1. We 
decompose any permutation a E Sk according to the image of 1, say a{l) = m. Upon 
relabelling indices, the corresponding permutation a' is in Sk-i, and we may apply the 
recursion hypothesis to the summation over such a'. We have 

k 

^(a,(u)) = 5^ (-1)'^A,(u-^ 



m=l crESk 
a{l)=m 

^_l).(/c-i)/2 A ^!h!l^±l^u,U2---Uk-2A{n) 

16 



(4.18) 



where we have apphed the recursion hypothesis to Aq{ui, . . . ,Um-i,Um+i, ■ ■ ■ ,Uk) and 
reabsorbed the sign change by (—1)'^"^ into the ratio A(u) / Ylt^mi'^^ ~ ^^)- conclude, 
we stiU have to prove the foUowing sublemma: 

k k 

ipk{ui,...,Uk) = } ^ = Uk-i (4.19) 

m=l i=l 

for ah distinct complex numbers mi, . . . , Wfc. To prove the latter, let us first note that it 
is a rational fraction, symmetric in the w's. Viewed as a function of ui, it has possible 
poles ai U2,U3, . . . ,Uk and is bounded at infinity. By symmetry, it is sufficient to compute 
the residue at Ui ^ W2, for which only the two first terms in the summation contribute, 
leading to: 

Res.,^., ( 17 _ <1^^ - g-^^i rr <in2-q-'uA ^ ^ ^^ ^0) 

Y Wl - W2 fj-^ U2 - Ui Ml - U2 fj-^ U2 - Ui J 

Hence the function (/9fc is bounded and has no pole in ui, hence is independent of wi, but 
as it is symmetric, it is a constant, say Ck- To compute it, we take the limit ui oo, 
and find the recursion relation (^fc(wi, . . . ,Uk) = q^~^ + q~^ipk-i{u2, ■ ■ ■ ,'U/c), henceforth 
Ck = + (l~^Ck-i- Moreover, by direct inspection, we find Ci = 0, therefore the 
sequence C is entirely fixed, and coincides with that of the Chebyshev polynomials of the 
second kind, namely Ck — Uk-i = {q^ — q~^)/{q — q~^)- This completes the proof of 

im)- 

Applying (|4.17| ) to ( [4.16| ), and decomposing when necessary the permutation a ac- 
cording to the images of 1 and/or k, say cr(l) = j and ct(/c) = m, we arrive at 



(- 



A(T{u^,...,Uk)) ( ^ \ 



k k _i 



- [Zi+i - Zi)< Uk-2Uk-3 - Uk-3 } Jj Uk-3 } Qm 

■'J.J. u ^-^ J. J. ^, _ 

1, j = l i=\ ■? ^ m=\ i=\ ^ ^ 

quj - q~'^Um Yj quj - q~^Ui quj - q~^Um 

/ , Jj9m I I 

i<j^<k ^^-"^ Li ^^-^^ 

(4.21) 

where we use the notations ( [4.15| ). Our last task is to prove that the r.h.s. of ( [4.21J ) 
vanishes identically (we denote it by B in the following). To this end, we view it as a 



m 
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rational fraction of the variable Zi+i, with possible poles at Zi+i Ug, £ = 1,2, . . . , k and 
at infinity. We first compute the residue at 2^+1 — > Ug. From the definition ( [4.15|) , we have 
KeSzi_^_J^^us{9s) = g's = ~{l ~ Q~^)^si and all other terms have a finite limit, henceforth: 

k k 2 



9 s t=\ t=i ^^-^'^ 




(4.22) 



J. J. 1l« 91 1 J- J- J- J- -))/, 



'k-3 

Uj — U£ 



Ui — Us — Uj — U£ 

e=i ^ ^ 1=1 j=i £=1 ^ ^ 



Noting that the last bracket involves only functions of the k—1 variables ui, . . . , Us-i, Us+i, ■ ■ ■ 
its vanishing is actually the consequence of the following lemma, valid for all p > 1 : 

flfe-tf^fl ''T'y +"--^ = ° (4.23) 



^=1 j = l £=l 



Uj — U£ 



Proof. Viewed as a function of Zi via the definition (|4.15|) , the l.h.s. of ( [4.23| ) (which we 
denote by D) is a rational fraction, with possible poles at Zi q^U£, i = 1,2, ... ,p and at 
infinity. Let us first compute the residue at Zi —>■ q^Us, for which the only contributions 
come from ReSz.^q2^^[fs) = {q - q~^)us = 

£^s £^s 

So D has no finite pole, and it is moreover bounded at infinity, with limit 

hm {D)=q^-qj2fl ^^^l^^ + U,., (4.25) 



j=l £=1 
£^P 



Uj — Ui 



Applying the above sublemma ( |4.19D , we may rewrite this into q^ — qllp-i + Up-2 = 0. We 
conclude that D = 0, and henceforth B has no finite pole in Zi+i. We must now examine 
possible residues at Zi+i oo. The leading behavior of B when Zi+i — > cxd is polynomial 
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of degree < 1. Noting that all lim^.^^^oo = Q ^, the coefficient of ZiJ^i of B in this limit 
reads: 

e=l j = l £=1 m=l 1=1 



(4.26) 

Applying the sublemma (f4.19| ) with /c — /c — 1 in the first line and k in the second, this 
simplifies into 



k k k 



-1. 



■^^ ■'J.J. 

£=1 j = l 1=1 ^ ^ 

(4.27) 

After noting that q~^Uk-2 — Uks = q^~^ and Uk-2 — <l~^Uk-i = we finally apply 

the lemma ( [1.23| ) with p = k, with the result 



BU^^, = q-^ - q'-''Uk-2 + q-^'Uk-s = ^-^g'^"' + Uks - qUk-2) = (4.28) 

Hence we conclude that B has no finite pole in Zi^i and is bounded at 2^+1 00, it is 
therefore independent of ^i+i. We now evaluate B at z^+i — 0, in which case all ge = q, 
and 



B\zi+i=o — -Zi 



^ q + q'-' n - {qUk-2 - t^.-s) ^ H . - U,.,{U,.2 - qU.-i} 



\ 



Uj — Up 

1=1 j=i e=i ^ * 



(4.29) 

where we have used again the sublemma ( [4.19|) . Again, we note that qUk-2 — Uk-3 = q^~^ 
and Uk-2 — qUk-i = —q'^ j and we apply the lemma ( [4.23| ) with p = k to get 

SU,+,=o = -z^{q - q''-^Uk-2 + q^Uk-s) = -z,q\q^-'' + U^-^ - q-^Uk-2) = (4.30) 

We may therefore conclude that 5 = identically, which implies that A(T) = = ^(5'), 
which in turn implies ( [4.10| ), as explained above. This completes the proof of Theorem 2. 
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4-3. Case of reflecting boundary conditions 

All that has been described above can for example apply to the case of periodic 
boundary conditions treated in P], avoiding the lengthy discussion found in this paper; in 
this case the function F is just 1. In the present case of reflecting boundaries, the solution 
to the qKZ equation must incorporate the new boundary conditions, cf Eqs. cj ), which 

correspond to a non-trivial choice of F, namely: 

"J n ^ 



X{l<l<m<ni'^rn " Wl){qWl - Q ^W^){q^ - WeWm) ni<^<m<n(9^ " ^^^^ 



(4.31) 



where we recall that the contours of integration encircle the Zi but not q "^Zi, nor q'^z'- 



By a computation of residues similar to what is performed in , it is easy to show that 
^ai,...,a„ is in fact a polynomial in the variables zi, . . . ,zn of degree 3n(n — 1). 

We now want to show that these identify with the components in the intermediate 
basis of the solution of the qKZ equation with reflecting boundary conditions discussed in 
Sect. 3.1. We note that both quantities satisfy the main equation ( |4.1C1| ) (exchange relation) 



of the gKZ system, as a consequence of Theorems 1 and 2 above. Furthermore, direct 
computation shows that when ai = i for i = 1, 2, ti, the only non-vanishing contribution 
to the integral ( [4.31| ) comes from the multiresidue at Wi Zi, for z = 1, 2, n. Cancelling 
all denominators, we are left with the polynomial: 

^i,2,...,n = JJ (qzi - q~^Zj){q'^ - ZiZj) JJ {q Zi - q~^ Zj){q'^ - ZiZj) (4.32) 

l<i<j<n n+l<i<j<2n 

which is nothing but \E'7ro ( |3.3| ), and indeed the change of basis implies that \E'i,2,...,n = ^ti-q- 
Moreover, introducing the lexicographic order on non-decreasing sequences (ai, . . . , a^), we 
note that the sequence (1, 2, n) is the smallest yielding a non-zero result for the integral 
( [1.31|) . Indeed, for any strictly smaller sequence, at least two residues will have to be taken 
at identical points say wg, Wm Zi, causing the result to vanish, due to the Vandermonde 
determinant of the w's in the numerator. 

The flnal step is to show that every component \E'a can be deduced from \I/i,2,...,n by 
use of Eq. ( |4.1U| ) . With respect to the lexicographic order on the non-decreasing sequences 
(ai, . . . , On), Eq. ( fl.lOp can be considered as a triangular linear system, allowing to compute 
^ai,...,a^_i,i,...,i,i+i,a^+fe,...,a„ terms of other \E'a with smaller index. It is easy to conclude 
fe-i 

from this that the \l/a are entirely determined by the non-zero component with smallest 
index (g). 
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4.4- Homogeneous limit 

Let us now evaluate the homogeneous limit of the components of by setting zi = 
Z2 = ■ ■ ■ = Z2n = I5 and by renormalizing it in such a way that = 1. This amounts 
to taking the integral formula ( [4.31J ), substituting Zi = 1 for all i, and dividing it out by 

^,,(1, 1, 1) = g3n(n-l)/2(^ _ ^-1^2n(n-l) ^ ^_ i ^ „(„_l)/2 _ 

Performing then the change of variables Wi = in the resulting expresion, we 

get: 



Jai,...,a„ (t) — T 



n(2n-l) 



n 

m=l 



dUr 



2nium 



<£<m<: 



\2n 



l<e<m<n 



Note that the integration contours encircle but not U2/U1 = (1 — t'^)/t. This integral 
formula is quite cumbersome and proves quite difficult to deal with. Instead, let us use 



the inversion property (|3.4| ) to express the homogeneous components as the reflected com- 
ponents ^'p(7r) at ^1 = 2:2 = ■ ■ ■ = Z2n = (T" ■ This simplifles drastically the pole structure 
of (|4.31|) , and finally, after division by ^'7rQ(l, 1, 1), leads to the following integrals: 



71 



dui 
2niu^„ 



Y\ (l-UlUm) 

l<e<m<n 



X 



Y\ [Um - W^) (1 + TUm + U^Um) [t + U^ + Um) 



l<i<m<n 



(4.33) 



We have the relation 



obtained by performing the change of variables Wi = 9 • 
Ja{iT) = -f^b(7r) for all n G LP2n, hence we may call the K^s the "arch closing" basis 
elements. 



The expression in brackets in formula (^4.33|) can be be thought of as a (polynomial) 
"quasi-generating function" for the Xb, that is the homogeneous limit of the ^'a/^i,2...,n- 
An important consequence of this formula is that the -ftTb are polynomials with integer 
coefficients in r. Since the change of basis from the intermediate basis to the basis of link 
patterns has entries that are also polynomials with integer coefficients in r, we conclude 
that the homogeneous '^■,^{t), normalized by \E' 



1, possess the same property, as 
had been conjectured in [0. In particular, at r = 1, where the \I/7r(l) are identified with 
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the ground state components of the loop model with reflecting boundaries normalized by 
\l/7r(,(l) = 1, we conclude that all these components are integers, as had been conjectured 
earlier in [^. 

Note that the same integrality argument applies equally well to the case of periodic 
boundary conditions treated in 0. 



5. Even case: proofs of various conjectures 

5.1. Solution at T ^ 

When T — »• 0, the integral is easily evaluated upon changing to variables Vm — 

Um/r and explicitly retaining only the leading terms when r ^ in each factor. This 
results in 

K,,,...,BAr)^r-'-^'^ L.J f[-^ n {v^-ve){l+v,+v^) (5.1) 



m=l l<e<m<n 



Upon identifying the product as the Vandermonde determinant A{v{l + v)) (using the 
notation A{z) = Y[i<e<m<ni^rn—Z£) for the Vandermonde determinant of z = {zi, . . . , Zn}, 



also equal to det(z^ )e,m=i,2,...,n): we may recast ( pTT| ) into a single determinant: 

(r) ~ T-'-^ ^ det ( I ^v-^-'^{l + 

l<e,m<n \ / Zmv 



~ r 



det C^-M (5.2) 

l<£,m<n \bi — m J 



The latter determinant is nothing but the number Nio{bi, • • • i ^n) of NILP in bijection 
with TSSCPP with fixed endpoints 6i,62,...,6n in their NILP formulation, while the 
power of r reads n{n — 1) — ~ 1) — Pi^)^ which is also the number of boxes in the 
box decomposition of the Dyck path associated to tt. 

This result may now be immediately translated into an estimate for the small r be- 
havior of the gKZ solution in the link pattern basis. Indeed, the change of basis ( [4.1| ) 
allows to identify "^-Kir) ~ i^rb(7r)(''') when r ^ 0. This is readily seen by writing 

v&^(r) = ^C-\r)^,«i^b(a)(r) (5.3) 

a 

and recalling that C~^(t) has all entries polynomial in r, and that it is lower triangular 
with respect to containment order of the Dyck paths associated to the link patterns, we 
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deduce that any tt such that C~^(r)7r,a is non-zero must be contained in a, hence have 
a strictly smaUer number of boxes if it is distinct from a. As K-^(^^^^{t^ behaves hke 
^n(n-i)-/3(a)^ we deduce that any contribution to the sum ( [5 .31 ) with a 7^ tt is subleading, 
as /3(a) < /3(7r). 

This completes the proof of the small r conjecture of Ref. 0] , in the case of even size, 
namely that 

^.{r) ~ r^^"^ iVio (&i(7r), . . . , 6n(7r)) (5.4) 
where 6i(7r) denote the positions of the arch closures of tt, counted from right to left. 

5.2. Solution at large r 

For large r, we obtain the leading contribution to (r) by changing variables to 

Vm = TUm in the integral formula ( [4.33| ), and retaining only the leading order in r within 
each factor in the integrand. This yields 

Kb.,...,bAr) ~ rEC'-i) I...lf[ -%Aiv) flil + v^f-' (5.5) 

By multilinearity of the Vandermonde determinant, this may be recast into 
i^.„...,.„(r) ~ rE(^-i) det ( f ^v^'^il + vf-') 

l<£,m<n \J ZmV / 

~rE(^-i) det f/~M (5.6) 

i<e,m<n \bi — m J 



To translate this into a result for '^t^{t), let us again consider the change of basis ([5.3|), 
and note that ~ 1) = ''^(^ ~ 1) ~ /3('^)- We wish to prove that '^'k{j) ~ K-^(^t^-^{t) 

at large r. 



Fig. 5: The strip decomposition of a typical Dyck path. To each ascending 
step (z — — 1) — > (i,/i) we associate a diagonal row oi h — 1 boxes as 
indicated, thus forming a strip of boxes. Strips associated to ascending steps 
exhaust all the boxes in the decomposition of the Dyck path. 
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The degree of the matrix element Ca^TriT) in r is given by the foUowing quantity. Define 
first /i(7r, a) as the sum over the arch openings of a of the total number of arches of n 
sitting above their position (an arch (z, 7r(z)) is said to sit above position j iS i < j < 7r(z)). 
The quantity h{n,a) is also the sum over the heights /ii(7r) in the Dyck path of tt (or 
equivalently the position occupied by the path on the integer half-line at time i), measured 
at the positions i of the points in the Dyck path of a reached by an ascending step (i.e. 
such that hi{a) = hi-i{a) + 1), namely: 

2n 

hi(a) = /ii_i(Q:) + l 

With this expression, it is easy to see that h{n, n) = f3{n) + n for any n G LP2n- Indeed, 
h{n,n) is the sum of heights of ends of ascending steps in the Dyck path of tt. As illustrated 
in Fig. 5, we may associate to each such ascending step (z — 1, /ii_i(7r)) (z, hi{7r)) with 
hi^Tf) = hi-i{7r) + l the diagonal strip of /ii(7r) — 1 boxes with centers at (z+£ — 1, hi{'rc)—i), 
I = 1,2, hi{iT) — 1, and this exhausts all boxes of tt. Such a "strip decomposition" was 
considered in Ref. We deduce that (3{it) = X]r;ii(7r)=/ii_i(7r)+i(^*(^) ~ ^{'^^ tv) —n 
as there are exactly n ascending and n descending steps in the Dyck path. Then, using 
the definition (|4.2| ) and the fact that the Chebyshev polynomials Um have degree m in r, 
we have 



rfa,7r = deg (^Ca,7r(r)) = h{Tr, a) - h{7v, tt) (5.^ 



Finally we need the following lemma: 

The quantity /7r(a) = /i(7r, a) + h{a,a), where a runs over the hnk patterns whose 
Dyck path is included in that of tt, reaches its maximum at a = tt only. 
Proof: let us show that fnioi) is a non-decreasing function with the size of a, namely the 
number of boxes in the decomposition of its Dyck path. Assume a' G LP2n differs from a by 
a single box say at positions i — l,i,i-(-lin the Dyck path formulation, with identical heights 
hj{a') = hj{a) + 2dj^i except at the position i. We see easily that h{a' , a') = h{a, a) + 1 
as the ascending step {i, 1) in a is replaced by {i — 1, i) in a' , and hi{a') = hi^i{a) + 1. 
Moreover, 

/i(7r, a') = h{7T, a) + hii^n) — /ii_|_i(7r) > /i(7r, a) — 1 (5.9) 

as hi{Ti) > hi+i(7r) — 1. We deduce that f-K^a') > /7r(a), hence that /7r(7r) > fjria) for all 
a whose Dyck path is included in that of n. Finally, tt is the unique maximum of /7r(a), 
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as is easily seen by removing a box from n say at position i and comparing heights in 
the Dyck paths. Indeed, we have say (/ii_i(7r), /ii(7r), /ii_(_i(7r)) = {m,m + l,m), and in 
the box-removed tc' we have (/ii_i(7r'), hi{'K'), hi+i{n')) = (m, m — 1, m). Hence h{n, n') = 
/i(7r, Tv) — {m + 1) + m and h{Tv' , tv') = h(Tv, tt) — (m + 1) + m, so that fn{^') = fn{T^) — 2. 

We deduce from the above lemma an upper bound on the degree da^n ( p-8| ): da^-n = 
/i(7r, a) — /i(7r, tt) < /i(7r, tt) — a) = (3{t{) — /3(a). We may now prove that the large 
r contribution to '^-k^t) is given by that to Kb(7r)(^)- This is done by induction on the 
(decreasing) number of boxes in a. Assume it is true for all tt with a strictly larger 
number of boxes than a. Then for all these tt, '^■nij) has degree equal to that of iirb(7i-)(T), 
namely n{n — 1) — [3{tx). Then the expression Ca,7r(r)\l/7r(T), for a C tt and n ^ a 
has degree da^-K + n{n — 1) — [3{tx) < n{n — 1) — /3(a) by the above inequality. But 
-^b(a)('7") = YliTT C'a,7r(T)^7r(T) has degree n(n — 1) —/3(a), hence it must be attained by the 
term tt = a in the sum, and we deduce that ^'a(T) ~ -f^b(Q:)(^) for large r. As the result 
holds trivially for the largest link pattern ttq (the only non-vanishing matrix element of C 
with this first entry is just C7rp,7rp(r) = 1), this completes the desired proof. 

5.3. Generalized sum rule 

The fundamental remark of Ref. ||] for the even case was that summing \l/a over a 
specific subset of "opening arch" basis elements, namely the set of arch openings a's such 
that ai = 2z — 1 — Ci, G {0, 1}, amounted to summing ^'^^ over the whole set LP2n- This 
was readily seen as a property of the change of basis Ca,n{T) (|4.2|). Due to a refiection 
symmetry property, an analogous statement may be derived for the "arch closing" basis 
elements. As a result, the above sum rule for X^Ti-gLPan obtained by summing the 

integrals Ki,^,b2,...,b„{T) over the 6^ = 2z - 1 - e^, with G {0, 1}. 

More precisely, let us consider 

K{t\r)^ E t^'^K,_,^,3-e,,...,2n-i-.Ar) (5.10) 

This quantity has many interesting specializations: (i) t = corresponding to the "maxi- 
mum" component i^i,3,...,2n-i(T) = ^7r,„„^(T), where iTmax connects the points {2i — l,2z) 
via little arches only (ii) t ^ oo corresponding to -ftri,2,4,...,2n-2(T) which is also equal to a 
single component, namely the link pattern tt'^^^ connecting (2i,2i -|- 1), and (l,2n); (iii) 
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t = 1 corresponding to the sum rule ^j^g^p^™ ~ -^(Ik)- From the definition ( [4 .331 ), it 
is readily computed as 



K{t\ 



T 



n 

m=l 



Y[ (1 - UiU, 



l<£<m<j 



l<i<m<n 



(5.11) 



We now make use of the following lemma, first conjectured in P] and then proved by 



Zeilberger in |2l]] . 



Yl (1 - UiUm)A Yl 



u 



2-2m 



W {l+TUm + 



UlUr 



l<e<m<r 



, m=l 



l<e<m<n 



<0 



(5.12) 



l<e<m<r 



where the subscript < means that we retain in the corresponding multiple Laurent series 
of the Um the terms with only non-positive powers, and the symbol A means as before 
antisymmetrization with respect to the permutations of variables: A{f){ui, . . . ,Un) = 



E.6S„(-l)"/(«a(i), ■ ■ ■ , u^(n)). Noting that A{u) = AiUl^i and that § A{f)g 



§ fA{g) , and applying the lemma to ( |5.11|) , we get 



K{t\T 



n 

m=l 
n 

n 

m=l 



2mUr 



n {r + u, + u^)A[\{u];;^{T + u^^r-^) 



dur 



l<£<m<n 
,1— m/'_ I „ — l\m — 1 



m=l 



2niUr 



-(1 + turr^h'^^ir + u-'r-'A{u{r + u)) 



= det 



l<£,m<n \ T 2'Kiu 



du 



(1 + tu)u'-'^{T + U-')'^-\t + U) 



(5.13) 



where in the last step we have used the multilinearity of the Vandermonde determinant to 
rewrite the whole multiple integral as a determinant of single integrals. This finally yields 



K{t\T) 

/Am(t|r) 



^ det {f£,m{t\T)) 



r 



2^+2m-3-2r 



£-1 

r-l 
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m — 1\ / m — 1 
r — m I \r + l— m 



(5.14) 



(a) (b) (c) 

Fig. 6: (a) The fundamental domain of a modified CSTCPP and its NILP 
description (blue and black paths) and (b) the correspondence (heights of the 
box piles) to triangular arrays of integers (c). Colors are related to weights: 
red means a weight of r, purple green x. With these conventions, (a) 
receives a weight x x (times a global factor of t'^~^ = t^), while (b) 
and (c) have the weight x^. 

In order to reconnect with the results of 0, it is convenient to rewrite K(t\T), by 
shifting all indices: i = j + l,m = i + l,r = s + l, and noting that the first row/column 
do not contribute to the determinant: 

K{t\T) = det 

l<i,j<n— 1 

For generic t, the polynomials K{t\T) correspond to a refined r, t-enumeration of 
CSTCPP^, which can be described as follows. In the NILP formulation of the CSTCPP^ 
the CSTCPP^ are described by paths made by two orientations of lozenges in a fundamen- 
tal domain of the CSTCPP^: these are the (colored) lozenges of types A, B on Fig. 6(a). 
In Ref. 0, these paths were viewed as pairs of paths sharing their arrival point, the paths 
on one side being one step longer than on the other side: the set of paths below the dashed 
line of Fig. 6(a) was identified as the NILP in bijection with TSSCPP (paths represented 
in black), while that above was viewed as an augmented one, with one more last step in 
each path (paths represented in blue, with the last step taking place in the strip just above 
the dashed line). 

At each step of the paths, lozenges of type A above the diagonal dashed line and B 
below (in red on the figure) are given a weight of r in K{t\T), except the last step of the 
longer paths (in the strip just above the dashed line), where a factor t is given to lozenges 
of type B; it is however more convenient to consider for this last step that the factor of r 
is replaced by 1/t (purple lozenges on the figure) up to global multiplication by 



E 



2i+2j-2s 



2i 



tr 



2j-s + l 



+ 



J 

2j-s 



(5.15) 
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Let us now discuss various specializations of K{t\T). 

As announced above, at t = 0, we find the maximal component 



^^^..(r) = i^(0|r)= det 

l<-i,j<.n — 1 



.2i+2j-2s 



J 

2i — s J \2j — s 



(5.16) 



This matches the expression conjectured in (conjecture 3, Eq. (5.1)). Likewise, at 
t — >^ oo, we find: 



(r) = lim —KUW) = det 

'"""^ ^ t^oo ^ ' ^ l<i,j<n-l 



y~^^2i+2j-25 + l 



(5.17) 



We also obtain at t = 1 the sum rule for the components of 



det 

l<i,j<n — 1 



E 



r 



2i+2j-2s 



2i 



2j-s + l 



+ 



J 

2j-s 



(5.18) 



which matches the conjectured expression found in Ref. M (Eq. (4.10)). 



Finally, one more identification is of some interest. Setting t = r ^ we find the 



generating function T„(x = r'^,1) of 



The latter has several intepretations. One of 



them is the following: consider triangular arrays of non-negative integers a^j, i,j > 1, 
i+j<n, with weakly decreasing rows and columns and such that an < n — i + 1 for all i. 
These arrays turn out to be in bijection with CSTCPP'^, see Fig. 6(b); to produce T„(x, 1) 
one gives a weight x to parts aij such that aij < j — I, see Fig. 6(c). Via the bijection this 
corresponds in terms of plane partitions to lozenges of type B that are below the diagonal. 
We now show that this is the same weight that is given to CSTCPP^ in K{t-^\t). When 
t = all red/purple lozenges get a weight of r. Call Uab the number of lozenges of type a 
in region b where a = A,B,C and 6 =t, | corresponds to above/below the diagonal dashed 
line. Then the weig ht for CSTCPP^ is t^^^t+^bx^ whereas the weight for triangular arrays 
is x"^^^ . Now the number of tiles of each orientation is fixed (independent of the choice of 
plane partition): in particular n^f + ri^i = 'n{n + l)/2 — 1. Furthermore the number of 
tiles of the first 2 types in each region is also fixed: ua] +nB] = n{n — l)/2. We conclude 
immediately that ua^ =2nBi+n — l, which allows us to identify the weights taking 

into account the prefactor t^~^ and the equality x = r^. 

Tn{x, 1) is also conjectured to be the x-enumeration of VSASM where a weight x is 
given to each pair of —Is (conj. 3.2 of [^). In the current framework there is no obvious 
explanation of this coincidence. 



28 



6. The odd case 
6.1. General solution 

The hnk patterns in size 2n + 1 are obtained bijectively from those in size one more 
2n + 2, by simply erasing the rightmost arch and suppressing the rightmost endpoint 
2n + 2, so that the opening point of the rightmost arch remains isolated and unmatched. 
The component of \1/ in size 2n + 1 can then be obtained from the corresponding one in 
size 2?7, + 2 by setting zn = 0. There is of course another bijection which consists on the 
contrary in erasing the leftmost arch and endpoint, so that the closing point of the leftmost 
arch remains unmatched. This time setting zi = oo allows to obtain in odd size 2n + 1 
from even size 2n + 2. Since the qKZ equation is not quite left-right symmetric these two 
routes lead to different expressions. 

Here we use the second route, having in mind that eventually we shall apply left-right 
symmetry as before to express everything in terms of closings instead of openings. We 
start with Eq. ( [4.31| ) at = 2n -|- 2 with ai = 1 and integrate over wi; we set 2:1 = cx) at 
this stage and obtain after shifting one step the indices of the a's, w's and z's: 

i^i^j^^jv m=l 



(6.1) 



that is, exactly the same expression as Eq. (|4.31| ), but in which now = 2n + l. It provides 



us with \E' in odd size in terms of the openings a^, with the convention that the closing of 
the former leftmost arch is unmatched (the opening, i.e. the leftmost point, being erased). 



We now use left-right symmetry in the form of the inversion property (|3.4|) , to obtain 
the component ^'^^ ^ defined in terms of the n closings hi counted from right to left, 
with the convention that the opening of the former rightmost arch is unmatched. We 
simply take Eq. ( |6.1| ) and substitute Oi with hi and Zi with q'^ / Zi. Note that hi < 2i — 1 
for closings to the right of the unmatched point and hi < 2i for closings to its left. We can 
now take the homogeneous limit Zi = 1 and obtain 



m=l 



Y\ (1 - U£Um] 
l<e.<m<n 



X 

i<e.<m< 



J]^ [Um - Ui) (1 + TUm + U£Um) {t + U£ + Um) 



(6.2) 
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6.2. Solutions at r ^ and large r 

Repeating the calculations of Sects. 5.1 and 5.2, and using the same changes of vari- 
ables, we arrive easily at the following formulas, respectively for r ^ 0: 

det h-') (6.3) 

^' l<£,m<n \b£ — m / 

and r —>■ oo: 

<,...,.„(r)~rE(^-i) det ( ^ J (6.4) 

l<.l,m<.n \ui — fn J 

Note that the determinants (|6.3|) vanish when some of the 6j's attains 2j, hence the 
correct small r behaviour is actually sub leading in those cases. The answer (|6.3D is accurate 



only if the unmatched point of the corresponding link pattern is at position 1 (in which 
case all 6^ < 2z — 1). In general, the correct answer depends on the position of this point, 
and requires a further expansion in powers of r, that we shall omit here. 

Note also that the size 2n + 1 large r solution (§j^) is related to the one in size 2n + 2 
(^) via: 

K^,...,,Jr) K^,^+,^_^+,^iT) (6.5) 

Indeed, as mentioned above, the odd size 2n + 1 link patterns tc' are obtained from those 
of size one more 2?! + 2 tt by removing the rightmost arch and leaving its beginning 
point unmatched, so arch closings indeed get shifted by one unit in this bijection, namely 
6i(7r) = 1, and bi{iT) = bi-i{iT') + 1 for z = 2, . . . , n + 1. This makes the connection between 
Eqs. O and (U) explicit. 

The results ( |6.3|) - (|6.4|) translate immediately into the corresponding identical small 
and large r estimates for \l/7r(T), as the change of basis from arch openings to link patterns 
is directly inherited from that of the even case 2n + 2. 

6.3. Generalized sum rule 

An analogous statement as that of Sect. 5.3 may be derived for the case of odd size 
271 + 1, in which we now have to sum over "arch closing" basis elements i^T^^ ^ with 
6i = 2z — e^, G {0, 1}. The generalized sum rule now reads: 

K'{t\T)^ J2 ^''^^^2-e„4-e„...,2n-e„(r) (6.6) 

e«e{o,i} 



and reproduces ^7reLP2n+i at t = 1. 
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Repeating exactly the same steps as those of Section 5.3, with now an extra factor of 
11^=1(1 + '^""1 + ^m)/^m, we easily arrive at: 



K'{t\T) = det 

l<£,7Ti<n 



du 



det (0£,^(t|r) + (t)i^rn+i{t\r)) 

l<£,m<n 



where 



du 



\t-i 



_ ^2£+2m-2r--4 



r — 



m — 1 \ / m — 1 
r + 1 — m/ \r + 2 — m 



It is straightforward to show that 



K'{t\T) = det {g,,^{t\T) + gt-i,^{t\T)) 

l<l,m<n 



, det igi,^it\T)) 

l<l,m<n 



where we have defined 



9£, 



m — 1\ I m — 1 
ri ]+t{ 

2m — r J \2m — 1 — r 



(6.7) 



(6.9) 



(6.10) 



The proof of ( |6.9| ) is obtained by performing a term-by-term identification of (f)e,m + <Pe,m+i 
with g£-i,m+ge,m, and then identifying the determinant of these entries with det((7£,m(t|r)) 
by column manipulations of the latter matrix. 





2 2 1 
1 




(a) 



(b) 



(c) 



Fig. 7: (a) The fundamental domain of a CSTCPP and its NILP description 
(blue and black paths) and (b) the correspondence (heights of the box piles) 
to triangular arrays of integers (c). The color code for weights here is: red=r, 
purple=i, and green=a:. 
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The NILP formulation of CSTCPP given in Ref. is similar to the even case, except 
this time the paths are viewed as two sets of paths of equal lengths (each in bijection with 
TSSCPP). The generating function K'(t\T) corresponds to the generating polynomial of 
these objects, with a weight r for the same type of tiles as in the even case, see Fig. 7, 
except in the last steps of one of the set of paths, where this weight is replaced by t. 

The expression in the second line of (|6.9|) for t = 1 matches exactly the conjectured 
result of Ref. (Eq. (4.2)), equal to the generating polynomial for weighted CSTCPP. 

Likewise, at t = one recovers the conjectured expression of Ref. (conjecture 4, 
Eq. (5.3)), and the limit t ^ oo yields the corresponding reflected link pattern (conjecture 
3, Eq. (5.2)). Note that the latter is identical to the sum rule K{t~^\t) in even size one 
less (2n) with t = r~^, due to simple identities relating K' and K. Therefore it is also 
equal to the generating function T^ij'^ , 1) of [l22|| . 

Finally at t = r we can identify K'{r\T) with the generating function T^ij"^, 0) deflned 
in p2|. Once again, this is no surprise since T„(x, 0) is the generating function for triangular 
arrays of non-negative integers a^j, i,j > 1, i + j < n, with weakly decreasing rows 
and columns and such that an < n — i for all i, which turn out to be in one-to-one 
correspondence with CSTCPP, see Fig. 7. As in the case of even size, one can check that 
the reflnements are the same, i.e. that the weight x given to lozenges of one of the three 
types that are below the diagonal is the same weight that is given to CSTCPP in X'(r|r) 
if one sets x = r"^. Note that once again simple identities show that K'{t\t) is equal to 
K{0\t), that is the component '^n^a^i'^) ^i^^ 2n. 

7. Conclusion 

In this paper, we have proved various conjectures regarding the minimal polynomial 
solution of the qKZ equation with reflecting boundaries with q generic, expressed in the link 
pattern basis. This was done by writing the solutions as multiple residue integrals, modulo 
a triangular change of basis. As both integrals and the change of basis are completely 
explicit, we therefore end up with an explicit formula for each component \E'7r(^i5 •••5 ^n)- 
We hope these expressions will help us address the full Razumov-Stroganov conjecture 
which gives a conjectural interpretation for each in the homogeneous case zi = ... = 
zat = 1 (and q = — e*'^/^), and hopefully come up with a more general combinatorial 
interpretation of the polynomials ^1/71(7") in the homogeneous case for generic q. Note that 
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we have now a numerical recipe for calculating the \l/7r(r), via the explicit inversion of the 
change of basis, and an explicit generation of the integrals. 

As stressed and proved in this paper, the above change of basis is independent of 
the details of the boundary conditions imposed in addition to the main exchange relation 

= (cj — r)\l/. These details are simply reflected by the insertion of some specific 
symmetric function F in the definition of the integrals. The techniques of the present 
paper may therefore presumably be adapted to include the other boundary conditions 
considered in [^, parametrized by the root systems of classical Lie algebras. 

As shown in generalizations of the Razumov-Stroganov sum rule have been ob- 
tained and proved in the case of the level 1 qKZ equation pertaining to higher rank {slk) 
algebras at specific values of the parameter q {q = —e'^ and q = —1). We believe the 
construction of the present paper may be generalized to these cases, and may lead to new 
combinatorial interpretations, such as generalizations of Plane Partitions. 

Another direction of future research is to revisit the model of crossing loops considered 
in p^j25| , p6| , which is based on the Brauer algebra, and obtain integral formulae in the 
same spirit as those of the present work. It would then be particularly interesting to 
understand their interrelation with the geometry of the Brauer loop scheme. 

Acknowledgments 

We acknowledge the support of European Marie Curie Research Training Networks 
"ENIGMA" MRT-CT-2004-5652, "ENRAGE" MRTN-CT-2004-005616, ESF program 
"MISGAM" and of ANR program "GIMP" ANR-05-BLAN-0029-01. We wish to thank 
J. de Gier and P. Pyatov for pointing out the appearance in our work of the generating 



functions of |]2^, as well as for sharing their current research with us. We also thank 
R. Kedem and C. Krattenthaler for useful discussions. 



33 



References 



[1] A. V. Razumov and Yu. G. Stroganov, Spin chain and combinatorics, J. Phys. A 34 
(2001) 3185-3190, |arXiv : cond-mat/001214li 

[2] M.T. Batchelor, J. de Gier and B. Nienhuis, The quantum symmetric XXZ chain 
at A = —1/2, alternating sign matrices and plane partitions, J. Phys. A 34 (2001) 
L265-L270, |arXiv : cond-mat 7010138^ . 

[3] A.V. Razumov and Yu.G. Stroganov, Combinatorial nature of ground state vector of 
0(1) loop model, Theor. Math. Phys. 138 (2004) 333-337; Teor. Mat. Fiz. 138 (2004) 
395-400, larXiv : math . C0/0104216| . 

[4] A.V. Razumov and Yu.G. Stroganov, 0(1) loop model with different boundary con- 
ditions and symmetry classes of alternating-sign matrices, Theor. Math. Phys. 142 
(2005) 237-243; Teor. Mat. Fiz. 142 (2005) 284-292, |arXiv : cond-mat/0108103| . 

[5] P. Pearce, V. Rittenberg and J. de Gier, Critical Q=l Potts model and Temperley- 
Lieb stochastic processes, |arXiv : cond-mat/010805l| ; P. Pearce, V. Rittenberg, J. de 
Gier and B. Nienhuis, Temperley-Lieb Stochastic Processes, J. Phys. A 35 (2002) 
L661-L668, |arXiv : math-ph/02090T7| . 

[6] P. Di Francesco, Totally Symmetric Self Complementary Plane Partitions and the 
quantum Knizhnik-Zamolodchikov equation: a conjecture, J. Stat. Mech. (2006) 
P09008, larXiv : cond-mat/0607499| . 

[7] P. Di Francesco, Open boundary quantum Knizhnik-Zamolodchikov equation and 
the weighted enumeration of plane partitions with symmetries, J. Stat. Mech. (2007) 
P01024, larXiv :math-ph/061101^ . 

[8] P. Di Francesco and P. Zinn- Justin, Quantum Knizhnik-Zamolodchikov equation, gen- 
eralized Razumov-Stroganov sum rules and extended Joseph polynomials, J. Phys. A 
38 (2006) L815-L822, |arXiv : math-ph/0508059| . 

[9] P. Zinn-Justin and P. Di Francesco, Quantum Knizhnik-Zamolodchikov, Totally Sym- 
metric Self-Complementary Plane Partitions and Alternating Sign Matrices, accepted 
for pubhcation in Theor. Math. Phys., |arXiv : math-ph/ 0703015] . 
[10] J. de Gier, Loops, matchings and alternating-sign matrices, Discr. Math. 298 (2005) 

365-388, larXiv :math . C0/0211285| . 
[11] W.H. Mills, D.P. Robbins and H. Rumsey, J. Combin. Theory Ser. A 34 (1983), 340- 
359. 

[12] P. Di Francesco, Inhomogenous loop models with open boundaries, J. Phys. A 38 

(2005) 6091-6120, |arXiv : math-ph/050403"^ . 
[13] E. Sklyanin, Boundary conditions for integrable quantum systems, J. Phys. A 21 

(1988) 2375-2389. 

[14] M. Jimbo, R. Kedem, H. Konno, T. Miwa and R. Weston, Difference Equations in Spin 
Chains with a Boundary, Nucl. Phys. B 448 (1995) 429-456, |arXiv : hep-th/9502060| . 



34 



[15] P. Di Francesco and P. Zinn- Justin, Around the RazumovStroganov conjecture: 
proof of a multi-parameter sum rule, Elec. J. Comb. 12 (1) (2005) R6, [arXiv :math- 



ph/0410061 



[16] P. Di Francesco, Boundary qKZ equation and generalized RazumovStroganov sum 
rules for open IRF models, J. Stat. Mech. (2005) P09004, |arXiv : math-ph/05090ll 



[17] F.A. Smirnov, A general formula for soliton form factors in the quantum sine-Gordon 
model, J. Phys. A 19 (1986) L575-L578. 

[18] LB. Frenkel and N. Reshetikhin, Quantum afGne Algebras and Holonomic Difference 
Equations, Commun. Math. Phys. 146 (1992) 1-60. 

[19] M. Jimbo and T. Miwa, Algebraic Analysis of Solvable Lattice Models (AMS, Provi- 
dence, 1995). 

[20] P. Di Francesco, O. Golinelli and E. Guitter, Meanders and the Temperley-Lieb alge- 
bra. Comm. Math. Phys. 186 (1997) 1-59, |arXiv : hep-th/9602025| . 

[21] D. Zeilberger, Proof of a Conjecture of Philippe Di Francesco and Paul Zinn-Justin re- 
lated to the qKZ equations and to Dave Robbins' Two Favorite Combinatorial Objects 



in pittp : / / wm . math . rutger s . edu/ ^zeilberg/ p j . html 



[22] D. Robbins, Symmetry Classes of Alternating Sign Matrices, |arXiv : math . CO/0008045 



[23] P. Di Francesco and P. Zinn-Justin, From Orbital Varieties to Alternating Sign Ma- 
trices, extended abstract for FPSAC'06, |arXiv:math-ph/0512047| . 

[24] P. Di Francesco and P. Zinn-Justin, Inhomogeneous model of crossing loops and 
multidegrees of some algebraic varieties, Commun. Math. Phys. 262 (2006) 459-487, 
|arXiv:math-ph/0412031| . 

[25] A. Knutson and P. Zinn-Justin, A scheme related to the Brauer loop model. Adv. in 
Math. 214 (2007) 40-77, |arXiv : math . AG/0503224 . 

[26] A. Knutson and P. Zinn-Justin, Brauer loop scheme and orbital varieties. 



35 



